Casimir Effect in the Presence of External Fields 
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O. 

CN . In this work the Casimir effect is studied for scalar fields in the presence of bound- 

> 

q \ aries and under the influence of arbitrary smooth potentials of compact support. In 



this setting, piston configurations are analyzed in which the piston is modeled by a 
potential. For these configurations, analytic results for the Casimir energy and force 
are obtained by employing the zeta function regularization method. Also, explicit 
numerical results for the Casimir force are provided for pistons modeled by a class 



of compactly supported potentials that are realizable as delta-sequences. These re- 
sults are then generalized to higher dimensional pistons by considering additional 
Kaluza-Klein dimensions. 

I. INTRODUCTION 



,8 



In 1948 it was shown by Casimir |8j that the presence of a boundary significantly modifies 
the vacuum structure of a quantum field. In particular, he computed the pressure between 
two perfectly conducting parallel plates due to the ground state of an electromagnetic field 
and determined that an attractive force exists between them. As a result, an enormous 
amount of literature has been produced which has been focused on the analysis of the 
effects that the geometry of the boundary and the boundary conditions have on the vacuum 



structure of quantum fields (see e.g 



14 



31 



33] and references therein) . In the vast 



majority of cases, however, research has been centered on the analysis of quantum systems 
endowed with ideal boundary conditions. The reason for such polarized interest lies in the 
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fact that the eigenvalues and the eigenf unctions of the relevant Laplace type operator with 
ideal boundary conditions are often explicitly known. This makes the analysis of quantum 
vacuum effects, as a result, less involved. 

Ideal boundary conditions, although of crucial theoretical importance for understanding 
important aspects of quantum vacuum effects, do not always provide an accurate description 
of physical systems. For instance, it is well known that a neutral metal plate does not 
constitute an impenetrable boundary for all the frequencies associated with a quantum 
field. For this reason different methods have been suggested in order to construct models 
that more closely characterize physical properties of real materials. 

One of the approaches used to more precisely describe quantum systems relies on the 
replacement of ideal boundary conditions with a boundary modeled by a potential. The ra- 
tionale behind this idea is that it is possible to describe physical characteristics of a boundary 
by choosing an appropriate potential function. This process reduces to substituting a ge- 
ometric boundary of the system with a non-dynamical external field. Along these ideas, 
the analysis of the Casimir effect in the setting of potentials modeling a boundary have 
appeared, for example, in [j], [igJ . 

It is worth mentioning that the importance of external potentials in quantum systems is 
not confined to the description of a non-deal boundary. In fact, they p.ay an important pie 



in describing background potentials resultin 
sphalerons 28J and electroweak Skyrmions [2 



rom c 



13 



assical solutions like monopoles |34l. |37|. 



19-22 



35 



36|]. 



Furthermore, the formalism needed to study the vacuum energy of scalar fields under 
;he influence of integrable potentials in unbounded Euclidean space has been developed in 
ll [ill 1^1 • It is the purpose of the current work to extend this analysis to include finite 
spatial volumes and additional Kaluza-Klein dimensions. 

The models considered in this paper can be formally reduced to the analysis of the Casimir 
effect on a one-dimensional finite interval with an arbitrary, sufficiently smooth, potential. 
Assuming two additional free dimensions this represents a piston configuration in flat space, 
where the piston is modeled by a potential with a sharply concentrated support. The 
arbitrary additional dimensions are represented by a smooth Riemannian manifold JV with 
or without boundary, which describes the cross-section of the piston. Piston configurations 
have received significant interest in recent years as they are often free of divergences which 
allows for an unambiguous prediction of forces Q. Prior research has concentrated on 
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infinitely thin pistons where different types of ideal boundary conditions were imposed; see, 



e.g 



15 



18 



23 



25 



29|. As previously mentioned, this work differs substantially from 
previous ones in that infinitely thin pistons are replaced by potentials modeling pistons of 
finite thickness. 

In this paper the zeta function regularization method is used, in which the Casimir energy 
and force are given by contour integrals that involve boundary values of unique solutions 
to initial value problems. This particular approach has been successfully applied to the 



evaluation of functional determinants 



26 



27] . for which explicit results can be obtained 



through purely analytical means and can be represented in closed form. In order to compute 
the Casimir energy and the corresponding force, however, numerical integration through 
suitable quadrature methods is necessary. 

The outline of the paper is as follows. In Section [IT] the spectral zeta function associated 
with a boundary value problem in one dimension is represented in terms of a contour integral 
and the formalism needed to perform the analytic continuation to a suitable domain in the 
complex plane is developed. The obtained meromorphic extension is exploited to analyze the 
Casimir energy for a massive scalar field on an interval under the influence of a background 
potential. The Casimir force on the piston is found for cases where the potential mimics a 
piston; see Equation (I2.32p . In Section ITO1 numerical results for the Casimir force are given 
for various potentials constructed from smooth, compactly supported functions, for which 
each can be realized as a delta- sequence. In Section HVl additional Kaluza-Klein dimensions 
are included into the analysis and the Casimir energy is determined. In this setting, results 
for the Casimir force are provided for cases in which the additional dimensions are chosen 
to be either 1 or I 2 . In the Conclusions, the results are summarized and additional studies 
along the lines of this work are outlined. 



II. MASSIVE SCALAR FIELD ON A ONE-DIMENSIONAL INTERVAL 

In this section a non-selfinteracting massive scalar field is considered on the interval 
I — [0, L] C M under the influence of a smooth background potential V(x). The one-particle 
energy eigenvalues \i of this system are determined by the differential equation 

d 2 



2 +y(x) + m 2 )<Mx) = A £ <Mx), (2-1) 



i 



augmented by boundary conditions which, for definiteness, are chosen to be of Dirichlet type 

0,(0) = ML) = . (2.2) 

In what follows, the mass parameter m is assumed large enough so that all the eigenvalues 
of the problem ( 12. ID and ( 12.21) are strictly positive. 

The spectral zeta function associated with the above boundary value problem is defined 

by 

oo 

Cr(*) = X>r> (2-3) 

£=0 

and, due to the asymptotic behavior of the eigenvalues, it is convergent in the semi-plane 
dt(s) > 1/2. In the framework of zeta function regularization, the Casimir energy Ec as of 
the system is encoded in the value of (i{s) at s — —1/2 jflll^lQ]. More precisely 

E Cas = \im^Ci(a-l) , (2.4) 



a^O 2 V 2 , 

where \x represents an arbitrary parameter with the dimension of a mass. Since the point 
s = —1/2 does not belong to the domain of convergence of (12. 4p . it is necessary to perform 
the analytic continuation of the series in Equation (12.31) to a neighborhood of s = —1/2. 
Here, the strategy employed consists in rewriting the series ( 12.31) as a contour integral using 



Cauchy's residue theorem |l0j and then utilizing the obtained integral representation as 
a starting point for the analytic continuation. This technique has been successfully used 
or the calculation of functional determinants in the setting described by Equation ( 12. ip in 



26 



23| • 



Following the ideas developed in [17|, |26|, |27| , instead of the boundary value problem (12. ip 



and ( 12.21) . the following equivalent initial value problem is to be considered, 
d 2 \ 

■— + V{x)+m 2 )u u {x) = v 2 u v [x) , M0) = 0, <(0) = 1 , (2.5) 

where v G (D. The eigenvalues A, of the original eigenvalue problem (12. ip with Dirichlet 
boundary conditions ( 12. 2ft are then determined as solutions to the transcendental equation 

u v {L) = . (2.6) 

Let us point out that the solutions to Equation (12. 5p are uniquely determined and define an 
analytic function of v. 
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?or dt(s) > 1/2, the zeta function ( 12. 3p can be represented in terms of a contour integral 

as rij 

(i(s) = ^-j dv(is 2 + m 2 y s ^u u (L) , (2.7) 

7 

where the contour 7 encloses all solutions to Equation (12. 6p . assumed to be on the positive 
real axis, in the counterclockwise direction. A deformation of the integration contour to the 
imaginary axis in ( 12. 7p leads to the following expression 

00 

C/00 = / dk(k 2 -m 2 )- s ^u ik (L) , (2.8) 

which is now valid for 1/2 < 3?(s) < 1. The analytic continuation of (12. 8ft to the region 
3^( s ) < 1/2 is obtained by adding and subtracting the large-fc asymptotic behavior of the 



solution u ik (L) [17 ]. 



The needed asymptotic behavior can be determined by applying a standard WKB tech- 
nique to the unique solution of the initial value problem 

d 2 



+ V(x) + k 2 u lk (x) = , u lk (0) = , i4(0) = 1 . (2.9) 



dx 1 

Although in the process of analytic continuation one only needs to be concerned with the 
exponentially growing part for large k, at this stage it is important to take into account both 
the exponentially growing and the exponentially decaying contributions in order to be able 
to correctly impose the initial condition in ( I2.9p . Following, e.g., 0, 30], it is convenient to 
introduce the auxiliary function 

S{x,k) = d x hiij> k {x) , (2.10) 

where ip k (x) satisfies 

d 2 \ 

+ V{x) + k 2 U fe (x) = 0. (2.11) 



dx 2 

By using (I2.10p in (12. lip , it is not very difficult to show that S(x, k) satisfies the differential 
equation 

S'(x,k) = k 2 + V{x) -S 2 (x,k) , (2.12) 
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where, here and in the rest of this work, the prime indicates differentiation with respect to 
the variable x. The asymptotic expansion of S(x, k) for large k can be written in the form 

oo 

S( x , k) = k-'Siix) , (2.13) 

i=-l 

where the asymptotic orders Si(x) are recursively determined by 

S^(x) = ±1 , S (x) = , Sx(x) = ±^P- , (2.14) 



S i+l (x) = T \ {s'lix) +J2S j (x)S i _ j (x)j 



The two different signs in (12.141) correspond to the exponentially growing and decaying 
solutions ipk{x) of (12.111) . Let S ± (x, k) denote the solutions of (12.121) corresponding to the 
different signs, the associated solutions of ( 12. lip have the form 



The original function of interest, namely Ui k (x), is obtained as a linear combination 

u ik (x) = A+exp | J dt S + (t, A;) | + A~ exp | J dt S~ (t,fc)| , (2.16) 

where the arbitrary coefficients A + and A~ can be found by imposing the initial condition 
in (12.91) and they read 

A+ = ~ A ~' A+ = s+(o,k)-s-(o,k) ' (2 ' 17) 

By using the result ( I2.17P in the expression (I2.16P the large- A; behavior of u^L) can be 
found to be 

Uik(L) = ^ /n u 1 _ fn ^ exp < J dt S + (t, k) > (l + E(k)) , (2.18) 



S+(0,k)-S-(Q,k) 



o 



where E(k) denotes exponentially decreasing terms as k — > oo. For the relevant quantity in 
the integral (F2T8|) . 

L 

\nu lk {L) = -\n(S + (0,k) - S~(0,k)) + J dt S + (t,k) 

o 

oo 

= -ln(2fc) + kL + ^2djk~ j , (2.19) 

3=0 
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where exponentially small terms have been omitted and the coefficients dj are defined as 

L 



d 2 j+i = I dx ££. + ] (x 



(2.20) 



and 



with Qj(0) defined through the cumulant expansion 



In 



1+ £St-,<0) 



fc=i 



£ 

i=l 



fii(0) 



y 2; 



For completeness, the first six coefficients dj are explicitly given by 

L 

d = , d 1= l -Jdt V(t) , d 2 = ~[V(L) + V{0)} 



(2.21] 



(2.22) 



V(L)-no)i-i 



dt V 2 {t) 



d 4 = -l[\/"(L) + \/"(0)] + i[y 2 (L)-V 2 



^ (3) (£) - V®(0)] - ^[V(L)V'(L) - V(0)V 



1 

16 



(2.23) 



dt V 3 (t) 



1 

32 



dt y(t)y"(t) . 



It is clear that by using (I2.14p and the definition (12. 20 p . an arbitrary number of coefficients 
can be determined by using an algebraic computer program. 

By adding and subtracting from the integral representation (12. 8p the leading N + 1 terms 
in the asymptotic expansion (12.241) . the zeta function can be represented as a sum of two 
terms 



(as). 



(2.24) 



where 11 



£(/)( s ) = ^11 1 dk {k 2 - ™ 2 )~ s 4; ( lnUifc(L) - kL + ln(2fc) - J^djArA > ( 2 - 25 ) 



s 



and 

N 



III 



c M (s) = s _^l J dk (A .2 _ m 2 )-A I kL - \n(2k) + J^d^ I . (2.26) 



3=0 

By construction, the function Q (s) is analytic in the region 9ft(s) > — (N + l)/2 and the 
meromorphic structure of (l as \s) is made manifest once the ^-integration is performed. 
More explicitly, 

where it is clear now that £r ( s ) represents a meromorphic function in the entire complex 
plane possessing only simple poles. The expression ( I2.24p together with (I2.25P and (I2.26P 
represents the desired analytic continuation of the spectral zeta function (12. 7p . 

For the purpose of computing the Casimir energy it is sufficient to choose N = 1 in the 
above expressions for Cf^( s ) an d (i as \s)- In this case, (^/\s) is analytic for 3ft(s) > — 1 and 
one can simply set s = —1/2 in (12. 25ft to obtain 

oo 

CpC-1/2) = ~Jdk (k 2 - ™ 2 ) 1/2 -^ jln ^(L) - kL + ln(2k) - ^ J , (2.28) 

m 

while, in the neighborhood of s = —1/2, (I2.27P gives 

C <">(-l/2 + a) = i?*±^ 
a 47r 

+ — r-2mvr + Lm 2 (ln4-l) + 4rf 1 (ln2-l)-2(2rf 1 + Lm 2 )lnml + 0(a) . (2.29) 

The explicit form of the Casimir energy for this system easily follows by substituting ( I2.28P 
and (I2.29P in the following expression 

E Cas = (-0 + \ + hi//) ResC/ (-0 + 0(a) , (2.30) 

obtained by expanding ( 12 Ah about a = 0. In the above formula and throughout the rest 
of this paper Res denotes the residue of the function and FP its finite part. It is evident 
from (I2.29p and (I2.30j) that the Casimir energy is, in general, not well defined because of 
the presence of the term Res(/(— 1/2). In order to overcome this problem, the system is 
interpreted in terms of a piston configuration where the piston itself is modeled by the 
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potential V(x). For this purpose, the potential V(x) is assumed to have compact support 
within the interval [0, L\. More precisely, it is assumed that V(x) does not vanish for 
x G [a — e, a + e] C [0,L]. According to this description, the point x = a represents the 
position of the piston. 

The asymptotic terms dj are expressed either in terms of boundary values of V(x) and 
its derivatives or as an integral of V(x) and its derivatives, therefore they are independent 
of the position a. It follows from the previous remarks that the Casimir force on the piston, 
defined in terms of the Casimir energy as 

d 

Fcas = —-Q^Ecas , (2-31) 

is a well defined quantity since ResCr(— 1/2) is, in this setting, independent of a. Hence, the 
explicit expression for the force is 

oo 

F Cas (a) = ^jdk (fc 2 -m 2 ) 1/2 ^^lnu ik (L) . (2.32) 

m 

It is worth pointing out that according to the above formula the magnitude and direction 
of the force is encoded in the boundary values of the solution to an initial value problem 
associated with an ordinary differential equation and no additional information is necessary. 

Despite the simplicity of (I2.32p information about the behavior of Fc as as a function of 
a can only be obtained through numerical integration techniques since the solution of (12. 9p 
is not explicitly known for an arbitrary V(x). In the following section the analysis of Fc as 
is provided for different types of potentials constructed from smooth, compactly supported 
functions. 

III. EXAMPLES: GAUSSIAN POTENTIALS 

It is clear, from (I2.32p . that in order to extract information about the Casimir force Fc as 
the evaluation of Uik(L) is necessary. An immediate numerical concern is that solutions to 
the differential expression (12 .9p for large k contain an exponentially increasing term of the 
type e kx . For this reason, to ensure accuracy in the numerical evaluation, discretization sizes 
must be chosen to be sufficiently small. This restriction is computationally costly but can 
be circumvented with relative ease in the following manner. Let Uik(x) = e kx ipik(x) in (12.91) . 
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The newly introduced function ipik(x) satisfies the initial value problem 
d 2 . d 



2k— + V(x))cp ik (x) = 0, <pik(0) = 0, ^(0) = 1. (3.1) 
The Casimir force in ( 12 .32 j) can, therefore, be written as 

oo 

Fcas{a) = ±-Jdk(k 2 - m 2 ) 1/2 |-A [e kx \n^ k (L)} , (3.2) 

m 

and since the exponentially growing term does not depend on a, this simplifies to 

oo 

Fcas{a) = ±- [dk(k 2 - m 2 ) 1/2 |- ^ln^ fc (L) . (3.3) 
2n J oa ok 

m 

The expression (13. 3p is now suitable for a numerical evaluation since the exponentially 
growing term has been dealt with analytically and, as a result, stringent tolerances on the 
discretization sizes have been alleviated. 

In the next subsections, the results for the Casimir force on Gaussian potentials centered 
at a and with support of extension e are presented. More specifically, background potentials 
of the form are considered, 



V(x) 



cxp 



-(x-a) 2 
^-(x-a) 2 



for \x — a\ < e, 



(3.4) 



otherwise 



For simplicity, in the following examples m = 0, L = 1, and e = 10 -4 . 

In order to obtain the Casimir force on the piston as a function of the position a from 
the expression ( 13.31) an adaptive second order Runge-Kutta method is used to integrate 
the differential equation ( 13. ip . with a potential of the form ( 13. 4p . This approximates the 
value of (fik{x) at x = L. For potentials of the type f!3.4[) this method can be shown to 
be convergent. More importantly, the discretization size can be determined such that the 
error in the approximation is within a user-specified tolerance. Upon successful calculation 
of (fik(L), standard second order centered differences are used to approximate the necessary 
derivatives in the integral formulation for the Casimir force. 

In the next step, the integral over the finite interval (m, M) is computed through the 
use of a symplectic integrator. The cutoff parameter, M, in the integral is determined in 
the following manner. Let I M be the calculated Casimir force up to M. Now, since the 
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integral is convergent, there exists a value of M such that the contribution of the integral 
beyond M is negligible. In effect, this can be determined as the first value of M for which 
\Im+i — Im\ < $ f° r some arbitrary 5. Here, 5 is chosen to be of the same order as the 
error obtained in the numerical integration of ( 13. lft . This efficient process is highly accurate, 
allows for large flexibility, and offers improved confidence in the results obtained. 

A. Positive and Negative Potential 

By setting 7] = 1 in ( 13 Ah one obtains a smooth, positive potential possessing a maximum 
at the point a = 1/2. This potential is depicted in the first graph of Figured! The Casimir 
force on the piston modeled by this potential is plotted in the second graph of Figure [TJ 
The Casimir force in this case is negative when the potential is close to the left boundary 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 01 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x a 



Figure 1: (a) The potential is shown over the domain of its support, (b) The Casimir force, F cas , 
is calculated using a standard potential of radius e = 10 -4 centered at a and 77 = 1. 

at x = 0, while it is positive when the potential is close to the right boundary at x — 1. In 
addition, the force vanishes at a = 1/2 as one would expect. This means that the positive 
potential considered above is always attracted to the closest boundary, making a = 1/2 a 
point of unstable equilibrium. 

By setting 77 = —1 in (13.41) one obtains the potential illustrated in the first graph of Figure 
EJ This potential is smooth, negative, and possesses a minimum at the point a = 1/2. The 
Casimir force associated with this potential has been plotted in the second graph of Figure 
[2j In this situation the behavior of the Casimir force is exactly opposite to the one found 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x a 



Figure 2: (a) The potential is shown over the domain of its support, (b) The Casimir force, F cas , 
is calculated using a standard potential of radius e = 10~ 4 centered at a and rj = — 1. 

for the positive potential. In particular, the negative potential is always repelled from the 
closest boundary and a = 1/2 is, in this point of stable equilibrium. 

B. Doubly-peaked Positive and Negative Potential 

The doubly-peaked positive potential is constructed by setting rj = 1/2 and by adding two 
potentials of the form (13.41) one with center at l/2 + e/2 and the other centered at 1/2 — e/2. 
The resulting potential, in the first plot of Figure El is smooth, positive, possessing two 
maxima and a minimum at a = 1/2. The Casimir force for this potential is reported in 
the second graph of Figure El The behavior of the force as function of the position of the 
potential is qualitatively similar to the case of the positive potential. In particular, the 
doubly-peaked positive potential is always attracted to the closest boundary. 

The doubly-peaked negative potential is constructed in the same way as the doubly- 
peaked positive potential by setting rj = —1/2. The resulting potential, characterized by 
two minima and one maximum at a = 1/2, is depicted in the first graph of Figure H] and the 
associated Casimir force is plotted in the second graph. Hence, the Casimir force in this case 
is analogous to the one found for the negative potential. This means that the doubly-peaked 
negative potential is always repelled from the closest boundary. 

It is not surprising that the behavior of the Casimir force on the doubly-peaked potentials 
matches the one found for the single potentials considered in the previous subsection. This 
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V(x)os 




Figure 3: (a) The potential is shown over the domain of its support, (b) The Casimir force, F c 
is calculated using a double potential of radius e = 10 -4 centered at a. 




0.1 0.2 0.3 0.4 0.5 06 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x a 



Figure 4: (a) The potential is shown over the domain of its support, (b) The Casimir force, F cas , 
is calculated using a double potential of radius e = 10~ 4 centered at a. 

should be expected since, effectively, the potentials have been evenly split while maintaining 
the total area under the curve equal to the single potentials. 

C. Mixed Potential 



The mixed potential, illustrated in the first plot of Figure [SI is obtained by adding two 
potentials of the form (13. 4j) : One with 77 = 1/2 and center at 1/2 — e/2 and the other with 
7] = —1/2, centered at 1/2 + e/2. The resulting Casimir force acting on this potential is 
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provided by the second plot in Figure [5j It is observed that the force on the potential 
is always negative, in contrast to the other cases considered. In other words, the mixed 
potential is repelled from the right boundary at x — 1 but attracted to the left boundary 
at x = 0. It is interesting to notice that the opposite signs of the potential have the net 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 5: (a) The potential is shown over the domain of its support, (b) The Casimir force, F c 
is calculated using a double potential of radius e = 10~ 4 centered at a. 



effect of eliminating the force throughout the majority of the interval except for the regions 
that are closer to the boundary. In the proximity of the right boundary the negative part 
of the potential becomes dominant and, therefore, the resulting Casimir force in that region 
resembles the one for negative potentials. Near the left boundary, instead, the positive part 
of the potential becomes dominant resulting in a force that behaves similarly to the one 
found in the case of positive potentials. 



IV. HIGHER DIMENSIONAL PISTONS 



In this section higher dimensional pistons modeled by potentials and constructed as 
product manifolds are studied. Let M be a D = (d + l)-dimensional manifold such that 
M = [0, L] x jV , where jV is a smooth d-dimensional Riemannian manifold representing 
the additional Kaluza-Klein dimensions. From the manifold M a piston configuration is 
obtained by modeling the piston itself with a smooth potential V(x) having support in the 
interior of the interval [0, L}. In this setting, the manifold JY represents the cross-section of 
the piston. 
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The dynamics of massless scalar fields propagating on M under the influence of the 
potential V(x) is described by the operator 

Sf=~-A^ + V{x), (4.1) 

where A j/ denotes the Laplacian operator on the manifold JV . The eigenvalue equation 

J?4 > (x,y) = \ 2 4 > (x,y) , (4.2) 
is separable and its solutions can, hence, be written as a product 

</>(x,y) = X(x)<p{y) , (4.3) 

where y denote the coordinates on the manifold JV and <p(y) represent the eigenfunctions 
of A jy satisfying 

- Ajr<p t (y) = rfmiy) ■ (4.4) 

By substituting the expression H4.3[) in the equation (]4.2p . and by setting A 2 = z/ 2 + rjj one 
can show that the functions X(x) are solutions to the equation 

dx 2 

The eigenvalues v appearing in (14.51) are uniquely determined once the boundary conditions 
for X(x) have been specified. As previously done in Section [III Dirichlet boundary conditions 
are imposed, namely, 

X u (0) = X V {L) = . (4.6) 
The spectral zeta function for the system under consideration is 

c(s) = E A ~ 2s = E^ 2 + ^ 2 )" s ' ( 4 - 7 ) 

which converges for dt(s) > (d + l)/2. 

It is clear, at this point, that since the eigenvalue equation ( 14. 2 p is separable, leading 
to (j4.4p and to (14.51) . the analysis of the Casimir energy for higher dimensional pistons 
and the associated Casimir force can be performed by using the methods described for the 
one-dimensional case. 



+ V(x) )X u {x) = v 2 X v {x) . (4.5) 
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Following the ideas developed in Section [TT1 the spectral zeta function (14. 7p is represented 
in terms of a contour integral as follows 

^) = 2^E/ d^ 2 + 4r^u,{L) , (4.8) 

1 7 

where u^(x) are the solutions to the initial value problem (cf. Section |H]) 

' + V(x)] u^x) = fi 2 u^(x) , u„(0) = , <(0) = 1 . (4.9) 



cfe 2 

Due to the presence of the manifold JV the analytically continued expression for the 
spectral zeta function ((s) will be written in terms of the spectral zeta function associated 
with A_yy, namely 

Ck(*) = 5>~ 2S ' 

which is well defined for 9ft(s) > d/2 and can be extended to a meromorphic function in the 
entire complex plane possessing only simple poles. The standard technique of adding and 
subtracting asymptotic terms is performed, in particular, 1^] 

t u \s) = ^Y.f dk ^ 2 - ^~ s Tk ln ^- (L) ~ kL+ H2k) - £ % p (411) 

t i l 3=0 J 

Here, C ( s ) can be proved to be well defined for (D — N — 2)/2 < 9fJ(s) < 1. By choosing 
N = D, C^K S ) becomes an analytic function in the neighborhood of s = —1/2 and can, 
therefore, be used for the computation of the Casimir energy with no further manipulations. 
Using the well-known meromorphic structure of the spectral zeta function associated with 
the Laplacian on JV {24], the relevant expression for £(s) at s — —1/2 reads 



C (/) ("l/2) 



( D \ 

- W ^ ( fc2 - ^ xl2 Tk ln uAL) - kL+ H2k) - £ I [ ' (413) 
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and 



<e»>(-l/2 + £ ) = I Aof-l)-lResCW ■ * C " (0) 



e I 47r v ' 2 2tt 




, 2 ^ r (§) " V 2 

+ ^ [CV(-l) + Cr(-l) (In 4- 1)] - iFP (.,-(-1/2) 
+ |[CV(0) + C,-(0)(.n4-2)] + |:^=^ 

+ ResCr ( J J - 7 - In 4 + 2 



FPC 



, 0(6) , (4.14) 

where \l/(x) represents the logarithmic derivative of the Euler gamma function. According 
to (12. 4p . the Casimir energy is obtained by adding the expressions (I4.13P and (14.141) and by 
multiplying by 1/2. It is worth noting that the above results are quite general and valid for 
an arbitrary smooth Riemannian manifold JV and for any smooth potential with compact 
support in [0, L]. More explicit results can only be found once the manifold JV and the 
potential V(x) are completely specified. Despite the lack of an explicit expression for the 
Casimir energy, one can still make some general remarks. From the results (I4.13P and (14. 14[) 
it is clear that the Casimir energy is, generally, divergent. However, since none of the terms 
in 114. 14ft depend on the variable a, the resulting Casimir force acting on the piston acquires 
contributions only from the finite part in (14.131) and is, hence, finite. 

This section will be concluded by discussing two specific examples. Consider the par- 
ticular cases for which jV = IR and JV = 1R 2 , respectively. In such cases the continuous 
spectrum resulting from the unrestricted dimensions can be integrated to obtain the follow- 
ing spectral zeta function densities 

CzxMW = £ fe%, (s - I) , (4.15) 



/ 47rr(s) V 2 
respectively 

WW = 4^Zi)C/(«-l). (4-16) 
Using the known meromorphic structure of (i{s) it is then easily verified that 

oo 
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and 

oo 

- wl™ o (-1) = is / *<* 2 ~ ^' 2 im • < 418 > 

m 

where, again, FP denotes the finite part of the Laurent series expansion. Numerical results 
for the above expressions can be obtained by following the same procedure described in 
the previous section. When jV is either K. or M 2 the Casimir force on pistons modeled by 
the type of potentials described in Section II III is qualitatively similar to the force found in 
one dimension for the same type of potentials. Since in this case the plots of the force as 
a function of the position a resemble the ones in Figures HHS^b), for brevity, they are not 
included here. 



V. CONCLUSIONS 



In this paper a method is developed to study the Casimir energy for massive scalar fields 
confined in finite volumes under the influence of smooth background potentials. The pivotal 
point of our method relies on rewriting the eigenvalues of the relevant boundary value problem 
as solutions to a transcendental equation related to an equivalent initial value problem; see 
Equations ( 12. 5ft and ( 12. 6ft . The starting point of our approach is the representation of the 
spectral zeta function associated with our models in terms of a complex integral. The 
analytic continuation of £(s) to a neighborhood of s = —1/2 was then achieved by adding 
and subtracting a suitable number of asymptotic terms from the integral representation. 
Although the form of the potential has been left unspecified, standard WKB techniques 
have allowed us to effectively compute the asymptotic expansions needed for the analytic 
continuation. The obtained analytic continuation of the spectral zeta function is, then, used 
in order to compute the Casimir force acting on a piston modeled by a smooth potential 
with compact support. In this case, it is found that while the Casimir energy is, in general, 
divergent the Casimir force on the piston is well defined. This can be immediately understood 
by noticing that the divergent terms in the energy are independent of the position of the 
piston. 

In the framework of pistons modeled by potentials it is found that the Casimir force can 
only be computed numerically once a potential has been specified. For this reason, several 
types of Gaussian potentials have been considered and plots of the Casimir force on the 
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piston as a function of the position a have been provided. We would like to stress, at this 
point, that in our examples Gaussian potentials have been chosen for illustrative purposes 
only. In fact, the results obtained in this work are completely general and are valid for any 
smooth potential with compact support. The numerical results obtained for the Casimir 
force in the various examples are consistent with our physical intuition, a fact that improves 
confidence in the presented analysis. 

This work can naturally be continued by considering the effect that different types of 
boundary conditions, such as Neuman, Robin, or mixed, have on the Casimir force. This 
analysis would follow the lines presented here without any major technical complications. In 
fact, the boundary conditions determine uniquely the constants A + and A~ in the asymptotic 
expansion of the functions in ( 12. 1 6f) . Different types of boundary conditions will lead to 
different expressions for A + and A~ but will keep the form of the asymptotic expansion 
( I2.24p unchanged. The analytic continuation of the spectral zeta function would, then, 
proceed in the same way as presented in this work. 

In addition, of particular interest is to study higher dimensional pistons modeled by 
potentials when the additional Kaluza-Klein manifold JV allows for the explicit knowledge 
of the eigenvalue associated with A^. Along these lines the authors are currently in the 
process of analyzing spherically symmetric and cylindrically symmetric configurations where 
angular momentum sums introduce additional technical and numerical complications. 
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